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Abstract 

The beauty production cross section for deep inelastic scattering events with at 
least one hard jet in the Breit frame together with a muon has been measured, 
for photon virtualities > 2 GeV^, with the ZEUS detector at HERA using 
integrated luminosity of 72pb~^. The total visible cross section is o"^f;(ep 
e jet /i X) = 40.9 ± 5.7 (stat.)]'^44(syst.) pb. The next-to-leading order QCD 
prediction lies about 2.5 standard deviations below the data. The differential 
cross sections are in general consistent with the NLO QCD predictions; however 
at low values of Q^, Bjorken x, and muon transverse momentum, and high values 
of jet transverse energy and muon pseudorapidity, the prediction is about two 
standard deviations below the data. 



The ZEUS Collaboration 

S. Chekanov, M. Derrick, J.H. Loizides^, S. Magill, S. Miglioranzi^, B. Musgrave, J. Repond, 
R. Yoshida 

Argonne National Laboratory, Argonne, Illinois 60439-4815, USA " 
M.C.K. Mattingly 

Andrews University, Berrien Springs, Michigan 49104-0380, USA 
N. Pavel 

Institut fiir Physik der Humboldt- Universitdt zu Berlin, Berlin, Germany 

P. Antonioli, G. Bari, M. Basile, L. Bellagamba, D. Boscherini, A. Bruni, G. Bruni, 
G. Cara Romeo, L. Cifarelli, F. Cindolo, A. Contin, M. Corradi, S. De Pasquale, P. Giusti, 
G. lacobucci, A. Margotti, A. Montanari, R. Nania, F. Palmonari, A. Pesci, L. Rinaldi, 
G. Sartorelli, A. Zichichi 

University and INFN Bologna, Bologna, Italy ^ 

G. Aghuzumtsyan, D. Bartsch, I. Brock, S. Goers, H. Hartmann, E. Hilger, P. Irrgang, 

H. -P. Jakob, O. Kind, U. Meyer, E. PauP, J. Rautenberg, R. Renner, A. Stifutkin, 
J. Tandler^, K.C. Voss, M. Wang 

Physikalisches Institut der Universitdt Bonn, Bonn, Germany 

D.S. Bailey^ N.H. Brook, J.E. Cole, G.P. Heath, T. Namsoo, S. Robins, M. Wing 
H.H. Wills Physics Laboratory, University of Bristol, Bristol, United Kingdom"^ 

M. Capua, A. Mastroberardino, M. Schioppa, G. Susinno 

Galabria University, Physics Department and INFN, Gosenza, Italy ^ 

J.Y. Kim, I.T. Lim, K.J. Ma, M.Y. Pac^ 

Ghonnam National University, Kwangju, South Korea ^ 

M. Helbich, Y. Ning, Z. Ren, W.B. Schmidke, F. SciuUi 

Nevis Laboratories, Golumbia University, Irvington on Hudson, New York 10027° 

J. Chwastowski, A. Eskreys, J. Figiel, A. Galas, K. Olkiewicz, P. Stopa, L. Zawiejski 

Institute of Nuclear Physics, Gracow, Poland * 

L. Adamczyk, T. Bold, I. Grabowska-Bold^, D. Kisielewska, A.M. Kowal, M. Kowal, J. 
Lukasik, M. Przybycieh, L. Suszycki, D. Szuba, J. Szuba^ 

Faculty of Physics and Nuclear Techniques, A GH- University of Science and Technology, 
Gracow, Poland 

A. Kotanski^, W. Slominski 

Department of Physics, Jagellonian University, Gracow, Poland 



I 



V. Adler, U. Behrens, I. Bloch, K. Borras, V. Chiochia, D. Dannheim^, G. Drews, 
J. Fourletova, U. Fricke, A. Geiser, P. Gottlicher^", O. Gutsche, T. Haas, W. Hain, 
S. Hillert^^, C. Horn, B. Kahle, U. Kotz, H. Kowalski, G. Kramberger, H. Labes, D. Lelas, 

H. Lim, B. Lohr, R. Mankel, I.-A. Melzer-Pellmann, C.N. Nguyen, D. Notz, A.E. Nuncio- 
Quiroz, A. Polini, A. Raval, U. Schneekloth, U. Stosslein, G. Wolf, C. Youngman, W. Zeuner 
Deutsches Elektronen- Synchrotron DESY, Hamburg, Germany 

S. Schlenstedt 

DESY Zeuthen, Zeuthen, Germany 

G. Barbagli, E. Gallo, C. Genta, P. G. Pelfer 
University and INFN, Florence, Italy ^ 

A. Bamberger, A. Benen, F. Karstens, D. Dobur, N.N. Vlasov^^ 

Fakultdt fiir Physik der Universitdt Freiburg i.Br., Freiburg i.Br., Germany ^ 

P.J. Bussey, A.T. Doyle, J. Ferrando, J. Hamilton, S. Hanlon, D.H. Saxon, I.O. Skillicorn 

Department of Physics and Astronomy, University of Glasgow, Glasgow, United King- 
dom 

I. Gialas 

Department of Engineering in Management and Finance, Univ. of Aegean, Greece 

T. Carli, T. Gosau, U. Holm, N. Krumnack, E. Lohrmann, M. Milite, H. Salehi, P. Schleper, 
T. Schorner-Sadenius, S. Stonjek^^, K. Wichmann, K. Wick, A. Ziegler, Ar. Ziegler 
Hamburg University, Institute of Exp. Physics, Hamburg, Germany 

C. Collins- Tootli^^, C. Fondas, R. Gongalo^^, K.R. Long, A.D. Tapper 

Imperial Gollege London, High Energy Nuclear Physics Group, London, United King- 
dom ™ 

P. Cloth, D. Filges 

Forschungszentrum Jiilich, Institut fiir Kernphysik, Jiilich, Germany 

M. Kataoka^^ K. Na gano, K. Tokushuku^^, S. Yamada, Y. Yamazaki 
Institute of Particle and Nuclear Studies, KEK, Tsukuba, Japan ^ 

A.N. Barakbaev, E.G. Boos, N.S. Pokrovskiy, B.O. Zhautykov 

Institute of Physics and Technology of Ministry of Education and Science of Kazakhstan, 
Almaty, Kazakhstan 

D. Son 

Kyungpook National University, Genter for High Energy Physics, Daegu, South Korea ^ 
J. de Favereau, K. Piotrzkowski 

Institut de Physique Nucleaire, Universite Gatholique de Louvain, Louvain-la-N euve, Bel- 
gium 



II 



F. Barreiro, C. Glasman^^, O. Gonzalez, L. Labarga, J. del Peso, E. Tassi, J. Terron, 
M. Zambrana 

Departamento de Fisica Teorica, Universidad Autonoma de Madrid, Madrid, Spain ^ 

M. Barbi, F. Corriveau, S. Gliga, J. Lainesse, S. Padhi, D.G. Stairs, R. Walsh 
Department of Physics, McGill University, Montreal, Quebec, Canada H3 A 2T8 

T. Tsurugai 

Meiji Gakuin University, Faculty of General Education, Yokohama, Japan ^ 

A. Antonov, P. Danilov, B.A. Dolgoshein, D. Gladkov, V. Sosnovtsev, S. Suchkov 
Moscow Engineering Physics Institute, Moscow, Russia ^ 

R.K. Dementiev, P.F. Ermolov, I.I. Katkov, L.A. Khein, LA. Korzhavina, V.A. Kuzmin, 

B. B. Levchenko, O.Yu. Lukina, A.S. Proskuryakov, L.M. Shcheglova, S.A. Zotkin 
Moscow State University, Institute of Nuclear Physics, Moscow, Russia ^ 

I. Abt, C. Biittner, A. Caldwell, X. Liu, J. Sutiak 
Max-Planck-Institut fiir Physik, Miinchen, Germany 

N. Coppola, G. Grigorescu, S. Grijpink, A. Keramidas, E. Koffeman, P. Kooijman, 
E. Maddox, A. Pellegrino, S. Schagen, H. Tiecke, M. Vazquez, L. Wiggers, E. de Wolf 
NIKHEF and University of Amsterdam, Amsterdam, Netherlands ^ 

N. Briimmer, B. Bylsma, L.S. Durkin, T.Y. Ling 

Physics Department, Ohio State University, Golumbus, Ohio 43210"' 

A.M. Cooper-Sarkar, A. Cottrell, R.C.E. Devenish, B. Foster, G. Grzelak, C. Gwenlan^^, 
T. Kohno, S. Patel, P.B. Straub, R. Walczak 

Department of Physics, University of Oxford, Oxford United Kingdom 

P. Bellan, A. Bertolin, R. Brugnera, R. Carlin, F. Dal Corso, S. Dusini, A. Garfagnini, 
S. Limentani, A. Longhin, A. Parenti, M. Posocco, L. Stanco, M. Turcato 
Dipartimento di Fisica dell' Universitd and INFN, Padova, Italy ^ 

E.A. Heaphy, F. Metlica, B.Y. Oh, J.J. Whitmore^^ 

Department of Physics, Pennsylvania State University, University Park, Pennsylvania 
16802° 

Y. Iga 

Polytechnic University, Sagamihara, Japan ^ 

G. D'Agostini, G. Marini, A. Nigro 

Dipartimento di Fisica, Universitd 'La Sapienza' and INFN, Rome, Italy ^ 

C. Cormack^o, J.C. Hart, N.A. McCubbin 

Rutherford Appleton Laboratory, Ghilton, Didcot, Oxon, United Kingdom ™ 



III 



C. Heusch 

University of California, Santa Cruz, California 95064, USA " 
I.H. Park 

Department of Physics, Ewha Womans University, Seoul, Korea 
H. Abramowicz, A. Gabareen, S. Kananov, A. Kreisel, A. Levy 

Raymond and Beverly Sackler Faculty of Exact Sciences, School of Physics, Tel-Aviv 
University, Tel- Aviv, Israel 

M. Kuze 

Department of Physics, Tokyo Institute of Technology, Tokyo, Japan ^ 

T. Fusayasu, S. Kagawa, T. Tawara, T. Yamashita 
Department of Physics, University of Tokyo, Tokyo, Japan ^ 

R. Hamatsu, T. Hirose^, M. Inuzuka, H. Kaji, S. Kitamura^^, K. Matsuzawa 
Tokyo Metropolitan University, Department of Physics, Tokyo, Japan ^ 

M. Costa, M.I. Ferrero, V. Monaco, R. Sacchi, A. Solano 
Universita di Torino and INFN, Torino, Italy ^ 

M. Arneodo, M. Ruspa 

Universita del Piemonte Orientale, Novara, and INFN, Torino, Italy ^ 
T. Koop, J.F. Martin, A. Mirea 

Department of Physics, University of Toronto, Toronto, Ontario, Canada M5S lAV^ 

J.M. Butterworth22, R. Hall- Wilton, T.W. Jones, M.S. Lightwood, M.R. Sutton^ C. Targett- 
Adams 

Physics and Astronomy Department, University College London, London, United King- 
dom ™ 

J. Ciborowski23, R. Ciesielski^^ , p. Luzniak^^, R.J. Nowak, J.M. Pawlak, J. Sztuk^^, 
T. Tymieniecka, A. Ukleja, J. Ukleja^^, A.F. Zarnecki 

Warsaw University, Institute of Experimental Physics, Warsaw, Poland ^ 

M. Adamus, P. Plucinski 

Institute for Nuclear Studies, Warsaw, Poland 

Y. Eisenberg, D. Hochman, U. Karshon M. Riveline 

Department of Particle Physics, Weizmann Institute, Rehovot, Israel 

A. Everett, L.K. Gladilin^^, D. Kgira, S. Lammers, L. Li, D.D. Reeder, M. Rosin, P. Ryan, 
A.A. Savin, W.H. Smith 

Department of Physics, University of Wisconsin, Madison, Wisconsin 53706, USA 



IV 



S. Dhawan 

Department of Physics, Yale University, New Haven, Connecticut 06520-8121, USA " 

S. Bhadra, CD. Catterall, S. Fourletov, G. Hartner, S. Menary, M. Scares, J. Standage 
Department of Physics, York University, Ontario, Canada M3J 1P3 " 



V 



^ also affiliated with University College London, UK 

^ retired 

^ self-employed 

^ PPARC Advanced fellow 

^ now at Dongshin University, Naju, South Korea 

^ partly supported by Polish Ministry of Scientific Research and Information Technology, 
grant no. 2P03B 12225 

partly supported by Polish Ministry of Scientific Research and Information Technology, 
grant no.2P03B 12625 

* supported by the Polish State Committee for Scientific Research, grant no. 2 P03B 
09322 

^ now at Columbia University, N.Y., USA 
1° now at DESY group FEB 

now at University of Oxford, UK 

partly supported by Moscow State University, Russia 

now at the Department of Physics and Astronomy, University of Glasgow, UK 

now at Royal Holoway University of London, UK 

also at Nara Women's University, Nara, Japan 

also at University of Tokyo, Japan 

Ramon y Cajal Fellow 

PPARC Postdoctoral Research Fellow 

on leave of absence at The National Science Foundation, Arlington, VA, USA 
now at Queen Mary College, University of London, UK 

present address: Tokyo Metropolitan University of Health Sciences, Tokyo 116-8551, 
Japan 

also at University of Hamburg, Alexander von Humboldt Fellow 
also at Lodz University, Poland 

supported by the Polish State Committee for Scientific Research, grant no. 2P03B 
07222 

Lodz University, Poland 

Lodz University, Poland, supported by the KBN grant 2P03B12925 
2^ supported by the KBN grant 2P03B12725 

on leave from Moscow State University, Russia, partly supported by the Weizmann 
Institute via the U.S. -Israel Binational Science Foundation 



VI 



supported by the Natural Sciences and Engineering Research Council of 
Canada (NSERC) 

supported by the German Federal Ministry for Education and Research 
(BMBF), under contract numbers HZIGUA 2, HZIGUB 0, HZIPDA 5, 
HZIVFA 5 

supported in part by the MINERVA Gesellschaft fiir Forschung GmbH, the 
Israel Science Foundation (grant no. 293/02-11.2), the U.S.-Israel Binational 
Science Foundation and the Benozyio Center for High Energy Physics 
supported by the German-Israeli Foundation and the Israel Science Foundation 
supported by the Italian National Institute for Nuclear Physics (INFN) 
supported by the Japanese Ministry of Education, Culture, Sports, Science 
and Technology (MEXT) and its grants for Scientific Research 
supported by the Korean Ministry of Education and Korea Science and Engi- 
neering Foundation 

supported by the Netherlands Foundation for Research on Matter (FOM) 
supported by the Polish State Committee for Scientific Research, grant no. 
620/E-77/SPB/DESY/P-03/DZ 117/2003-2005 

partially supported by the German Federal Ministry for Education and Re- 
search (BMBF) 

supported by RF President grant N 1685.2003.2 for the leading scientific 
schools and by the Russian Ministry of Industry, Science and Technology 
through its grant for Scientific Research on High Energy Physics 
supported by the Spanish Ministry of Education and Science through funds 
provided by CICYT 

supported by the Particle Physics and Astronomy Research Council, UK 
supported by the US Department of Energy 
supported by the US National Science Foundation 

supported by the Polish Ministry of Scientific Research and Information Tech- 
nology, grant no. 112/E-356/SPUB/DESY/P-03/DZ 116/2003-2005 
supported by the Polish State Committee for Scientific Research, grant no. 
115/E-343/SPUB-M/DESY/P-03/DZ 121/2001-2002, 2 P03B 07022 



VII 



1 Introduction 



Deep inelastic scattering (DIS) offers a unique opportunity to study the production mech- 
anism of bottom (6) quarks via the strong interaction in a clean environment where a 
point-like projectile, a photon with a virtuality Q^, collides with a proton. Due to the 
large centre-of-mass energy, bb pairs are copiously produced at the electron-proton col- 
lider HERA. The large 6-quark mass provides a hard scale, making perturbative Quantum 
Chromodynamics (QCD) applicable. However, a hard scale can also be given by the trans- 
verse jet energy and by Q. The presence of two or more scales can lead to large logarithms 
in the calculation which can possibly spoil the convergence of the perturbative expansion. 
Precise differential cross-section measurements are therefore needed to test the theoretical 
understanding of 6-quark production in strong interactions. 

The cross sections for 6-quark production in strong interactions have been measured in 
proton-antiproton collisions at the SppS [1] and the Tevatron [2] and, more recently, in 
two-photon interactions at LEP [3] and in 7p interactions at HERA [4,5]. Some of the 
6-production cross sections are significantly above the QCD expectations calculated to 
next-to-leading order (NLO) in the strong coupling constant, a^. 

This paper reports the first measurement of 6-quark production in DIS at HERA, in the 
reaction with at least one hard jet in the Breit frame [6] and a muon, from a b decay, in 
the final state: 

ep^ebbX-^e+jet + fi + X. 

In the Breit frame, defined by 7 + 2xP = 0, where 7 is the momentum of the exchanged 
photon, X is the Bjorken scaling variable and P is the proton momentum, a space-like 
photon and a proton collide head-on. In this frame, any final-state particle with a high 
transverse momentum is produced by a hard QCD interaction. 

In this paper, a measurement of the visible cross section, cr^^, is presented, as well as 
several differential cross sections. The measured cross sections are compared to Monte 
Carlo (MC) models which use leading order (LO) matrix elements, with the inclusion of 
initial- and final-state parton showers, as well as to NLO QCD calculations. All cross 
sections are measured in a kinematic region in which the scattered electron, the muon 
and the jet are well reconstructed in the ZEUS detector. 

2 Experimental conditions 

The data used in this measurement were collected during the 1999-2000 HERA running 
period, where a proton beam of 920 GeV collided with a positron or electron beam of 
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27.5 GeV, corresponding to a centre-of-mass energy of 318 GeV. The total integrated 
luminosity was (72.4 ± 1.6) pb~^. 

A detailed description of the ZEUS detector can be found elsewhere [7,8]. A brief outline 
of the components that are most relevant for this analysis is given below. The high- 
resolution uranium-scintillator calorimeter (GAL) [9] consists of three parts: the forward 
(FGAL), the barrel (BGAL) and the rear (RGAL) calorimeters. Each part is subdivided 
transversely into towers and longitudinally into one electromagnetic section (EMG) and 
either one (in RGAL) or two (in BGAL and FGAL) hadronic sections (HAG). The smallest 
subdivision of the calorimeter is called a cell. The GAL energy resolutions, as measured 
under test-beam conditions, are a{E)/E = 0.18/ a/-E (GeV) for electrons and a{E)/E = 
0.35/v/^(GeV) for hadrons. 

Gharged particles are tracked in the central tracking detector (GTD) [10], which operates 
in a magnetic field of 1.43 T provided by a thin superconducting solenoid. The GTD 
consists of 72 cylindrical drift-chamber layers, organised in nine superlayers covering the 
polar-angle^ region 15° < 6 < 164°. The transverse-momentum resolution for full-length 
tracks can be parameterised as (t{pt)/pt = 0.0058^^ © 0.0065 © 0.0014/^7-, with px in 
GeV. 

The position of electrons^ scattered at small angles to the electron beam direction was 
measured using the small-angle rear tracking detector (SRTD) [11, 12]. The SRTD is 
attached to the front face of the RGAL and consists of two planes of scintillator strips, 
1 cm wide and 0.5 cm thick, arranged orthogonally. 

The muon system consists of tracking detectors (forward, barrel and rear muon chambers: 
FMUON [8], B/RMUON [13]), which are placed inside and outside a magnetised iron yoke 
surrounding the GAL and cover polar angles from 10° to 171°. The barrel and rear inner 
muon chambers cover polar angles from 34° to 135°. 

The luminosity was measured from the rate of the bremsstrahlung process ep — ^ e'-fp. The 
resulting small-angle energetic photons were measured by the luminosity monitor [14], a 
lead-scintillator calorimeter placed in the HERA tunnel at Z = —107 m. 



The ZEUS coordinate system is a right-handed Cartesian system, with the Z axis pointing in the 
proton beam direction, referred to as the "forward direction", and the X axis pointing left towards 
the centre of HERA. The coordinate origin is at the nominal interaction point. 
^ Hereafter "electron" refers both to electrons and positron. 
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3 Event Selection 



Events were selected online via a three-level trigger system [8, 15]. The trigger required a 
localised energy deposit in the EMC consistent with that of a scattered electron. At the 
third level, where a full event reconstruction is available, a muon was required, defined by 
a track in the CTD loosely matching a track segment in the inner part of the B /RMUON 
chambers. 

The scattered electron candidate was identified from the pattern of energy deposits in the 
CAL [16]. The energy (Ee) and polar angle {6e) of the electron are measured by combining 
the impact position at the calorimeter with the event vertex. The impact position is 
measured from the calorimeter cells associated with the electron candidate, but the CTD 
(6'e < 157°) and SRTD (6'e > 162°) detectors are used to improve the measurement 
whenever the electron trajectory lies within the respective regions of acceptance. 

The reconstruction of was based on the measurement of the scattered electron energy 
and polar angle [17]. The Bjorken scaling variables x and y were reconstructed using the 
E-method, which allows the determination of the estimator ?/s independently of initial 
state photon radiation by reconstructing the incident electron energy [18]. 

Events were selected [19] by requiring the presence of at least one muon in the final state 
and at least one jet in the Breit frame. The final sample was selected in four steps: 

1. Inclusive DIS event selection: 

• a well reconstructed scattered electron was required with energy greater than 
lOGeV, > 2 GeV^, yjB > 0.05 and < 0.7, where yjB is the y variable 
reconstructed using the Jacquet-Blondel method [20]; 

• for events with the scattered electron reconstructed within the SRTD acceptance 
the impact position of the electron was required to be outside a box defined by 
\Xe\ < 12 cm and \Ye\ < 6 cm. For events without SRTD information, a box cut 
on the face of the RCAL of |Xe| < 12 cm and \Yf;\ < 10 cm was used. This cut 
removed electron candidates near the inner edge of the RCAL beampipe hole; 

• to reduce the background from collisions of real photons with protons (photo- 
production), where the scattered electron escapes down the rear beampipe, the 
variable E — pz was required to be in the range 40 < E — pz < Q5 GeV. The 
variable E — pz was defined as the difference of the total energy and the longi- 
tudinal component of the total momentum, calculated using final-state objects, 
reconstructed from tracks and energy deposits in the calorimeter; 

• the reconstructed event vertex was required to lie within 50 cm of the nominal 
interaction point. 
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2. Muon finding: 

Muons were identified by requiring a track segment in both the inner and outer parts 
of the BMUON or RMUON chambers. The reconstructed muons were matched in 
space and momentum with a track found in the CTD, with a probability greater 
than 1%. This cut rejected the background from muons coming from K"^ and 
decays and from particles produced in hadronic showers in the CAL that may be 
misidentified as muons. In addition, cuts on the muon momentum, p'^, the muon 
transverse momentum, pj. and the muon pseudorapidity, rj^, were applied: 

• —0.9 < 7]^ < 1.3 and pf^ > 2 GeV corresponding to the BMUON region; 

• — 1.6 < ?7'* < —0.9 and > 2 GeV corresponding to the RMUON region. 

The reconstruction efficiency of the muon chambers was calculated separately for 
BMUON and RMUON using an independent data sample of di-muon events produced 
in photon-proton collisions [21]. This data sample consisted of elastic and quasi-elastic 
Bethe-Heitler events (77 and J/ip production. 

The data sample was selected from events triggered by the inner muon chambers. Two 
tracks, reconstructed in the CTD, with transverse momentum greater than 1 GeV and 
associated with energy deposits in the CAL consistent with a minimum-ionising parti- 
cle were required. One of the CTD tracks was required to point to the muon chamber 
that triggered the event, and the other was used to measure the muon efficiency, de- 
fined as the ratio of the number of tracks satisfying the muon matching requirement 
to the total number of tracks. The measured muon-reconstruction efficiencies are be- 
tween 20% and 40%, depending on the region of the muon chambers and on the muon 
transverse momentum. 

3. Jet finding: 

Hadronic final-state objects were boosted to the Breit frame and clustered into jets 
using the kx cluster algorithm (Ktclus) [22] in its longitudinally invariant inclu- 
sive mode [23]. The four-momenta of the hadronic final-state objects were calculated 
from the measured energy and angle, assuming the objects to be massless. The pt 
recombination scheme was used. Reconstructed muons were included in the cluster- 
ing procedure. Events were required to have at least one jet with transverse energy 
measured in the Breit frame, -E-Pjet*) above 6 GeV and within the detector acceptance, 
—2 < < 2.5, where rjj'^^ is the jet pseudorapidity in the laboratory frame. 

4. Muon-jet association: 

The muons in the sample were associated with the jet containing the corresponding 
hadronic final-state object using the Ktclus information. The associated jet was not 
necessarily the jet satisfying the jet requirements above. To ensure that the associated 
jet was well reconstructed, it was required to have -Erjet* > 4 GeV. 
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After these selection cuts, 941 events remained. 



4 Monte Carlo simulation and NLO QCD calcula- 
tions 

To correct the results for detector effects and to extract the fraction of events from b 
decays, two MC simulations were used: Rapgap 2.08/06 as default and Cascade 
1.00/09 for systematic checks. The predictions of the MC simulations were also compared 
to the final results. 

The program Rapgap 2.08/06 [24] is an event generator based on leading-order (LO) 
matrix elements, with higher-order QCD radiation simulated in the leading-logarithmic 
approximation using initial- and final-state parton showers based on the DGLAP equa- 
tions [25]. To estimate the background, samples with light and charm quarks in the final 
state were produced. The process in which a bb pair is produced in photon-gluon fusion 
was used to simulate the signal. The charm and 6-quark masses were set to l.SGeVand 

5 GeV, respectively. The CTEQ5L [26] parameterisation of the proton parton densities 
was used. Heavy-quark hadronisation was modelled by the Bowler fragmentation func- 
tion [27]. The rest of the hadronisation was simulated using the Lund string model [28] 
as implemented in Jetset 7.4 [29]. The Rapgap MC includes the LO electroweak 
corrections calculated using Heracles 4.6.1 [30]. 

The Cascade 1.00/09 MC [31] uses the 0{as) matrix elements, where the incoming 
partons can be off-shell. The parton evolution is based on the CCFM equations [32], 
which are derived from the principles of factorisation and colour coherence. The mass 
of the b quark was set to 4.75 GeV. 

The NLO QCD predictions were evaluated using the HvQDis program [33,34], which 
includes only point-like photon contributions. The fragmentation of a quark into a B 
hadron was modelled by the Kartvelishvili function [35]. The parameter a was set to 27.5, 
as obtained by an analysis [36] of e~^e~ data [37]. The semi-leptonic decay of B hadrons 
into muons was modelled using a parameterisation of the muon momentum spectrum 
extracted from Rapgap. This spectrum corresponds to a mixture of direct (b — > fi) and 
indirect {b ^ c fi) 5-hadron decays. Jets were reconstructed by running the inclusive 
kr algorithm, using the Pt recombination scheme, on the four-momentum of the two or 
three partons generated by the program. The 6-quark mass was set to = 4.75 GeV and 



the renormalisation and factorisation scales to yU = y -|- m^, where px.b is the mean 

transverse momentum of the b and b quarks. The CTEQ5F4 proton parton densities [26] 
were used. The sum of the branching ratios of direct and indirect decays of B hadrons 
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into muons was fixed to the Jetset 7.4 value of 0.22. 

Tlie NLO QCD predictions were multiplied by hadronisation corrections to compare them 
to the measured cross sections. The hadronisation corrections are defined as the ratio of 
the cross sections obtained by applying the jet finder to the four-momenta of all hadrons, 
assumed to be massless, and that from applying it to the four-momenta of all partons. 
They were evaluated using the Rapgap program; they lower the NLO QCD prediction 
by typically 10%. 

The uncertainty of the NLO prediction was estimated by varying the factorisation and 
renormalisation scales, /i, by a factor of 2 and the 6-quark mass, mb, between 4.5 and 
5.0 GeV and adding the respective contributions in quadrature. Additional uncertainties 
due to different scale choices and to different fragmentation functions are within the 
quoted uncertainties. More details of the NLO QCD calculation and of the determination 
of its uncertainties can be found elsewhere [33,34,38]. 

5 Extraction of the beauty fraction 

A significant background to the process under study is due to muons from in-flight decays 
of pions and kaons. Such decay muons are mostly characterised by low momenta and, 
therefore, partly rejected by the cuts p'^ > 2 GeV and pl^ > 2 GeV. In addition, the signal 
reconstructed in the muon chambers can be due to kaons or pions passing through the 
CAL. Muons can also originate from the semi-leptonic decay of charmed hadrons. These 
decays produce events topologically similar to those under study. 

Due to the large 6-quark mass, muons from semi-leptonic b decays usually have high 
values of the transverse momentum, p^^^ with respect to the axis of the closest jet. For 
muons from charm decays and in events induced by light quarks, the p^^' values are low. 
Therefore, the fraction of events from b decays in the data sample can be extracted on a 
statistical basis by fitting the relative contributions of the simulated bottom, charm and 
light-quark decays to the measured p^"' distribution. 

The extraction of the fraction of 6-quark decays relies on the correct simulation of the 
shape of the p^^' distribution for all processes. The simulation was checked with the data. 
For this purpose, an inclusive DIS data sample with at least one hard jet in the Breit frame 
was selected, without requiring a muon in the final state. For tracks passing the same 
selection criteria as required for the muon, the p^?' distribution was calculated. Figure CJi 
shows the comparison of the shape of the measured p^?^ distribution with the simulated 
light- and charm-quark contribution. The shape is reasonably well described. 

Figure [Dd shows the measured p^?^ distribution for muon candidates compared to the 
MC simulation. The MC simulation contains the background processes from light and 
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charm quarks and the contribution from b quarks. The distributions are peaked at low 
values, where the decays of hadrons containing charm and light quarks dominate. At 
higher values, the measured distribution falls less steeply than that expected for light- 
quark and charm contributions alone. To determine the 6-quark fraction in the data, the 
contributions from light-plus-charm flavours and beauty in the simulation were allowed 
to vary, and the best fit was extracted using a binned maximum-likelihood method. The 
measured fraction of events from b decays, is (30.2±4.1)%, where the error is statistical. 
The mixture with the fitted fractions describes the data well. 

Figures Ic-f show the comparison between the data and the MC simulation with respect to 
the momentum and the pseudorapidity of the muon, as well as the associated jet transverse 
energy in the Breit frame and the pseudorapidity of the associated jet measured in the 
laboratory frame. The MC simulation, with the different contributions weighted according 
to the fractions found using the fit procedure described above, reproduces the muon and 
jet kinematics well. 



6 Systematic uncertainties 

The systematic uncertainties on the measured cross sections were determined by changing 
the selection cuts or the analysis procedure in turn and repeating the extraction of the 
cross sections. The numbers given below refer to the total visible cross section, (7^5. For 
the differential distributions the systematic uncertainties were determined bin-by-bin and 
are included in the figures and in Table ^ The following systematic studies were carried 
out: 

• selection cuts and SRTD alignment: variation of the selection cuts by one standard 
deviation (including the electron energy, E - pz, E^'^-'^^, r]]^}' and SRTD box cut). 
This led to a systematic deviation of +9.1% and —6.1% with respect to the nominal 
value, where the biggest uncertainties were introduced by the widened rjjl^ cut and 
the increased E^^et cut. The relative alignment between the RCAL and the SRTD 
detector is known to a precision of ±1 mm [39]. The related systematic uncertainty 
was estimated by shifting the reconstructed SRTD hit position by ±2 mm in both 
coordinates and was +0.5% and —1.3%, respectively; 

• energy scale: the effect of the uncertainty in the absolute CAL energy scale of ±2% 
for hadrons and of ±1% for electrons was +3.3% and —0.3%; 

• extraction of b decays: the uncertainties related to the signal extraction were estimated 
by doubling and halving the charm contribution. This leads to a systematic uncer- 
tainty of +5.7% and —3.5%, respectively. The uncertainty obtained by reweighting 
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the light-plus-charm quark p^^' distribution with the one extracted from the data as 
described in Section El is within this uncertainty; 

• muon reconstruction efficiency: the effect of the uncertainty on the muon reconstruc- 
tion efficiency for the barrel and rear regions of the muon detectors was +8.9% and 
-7.8%; 

• model dependence of acceptance corrections: to evaluate the systematic uncertain- 
ties on the detector corrections, the results obtained with Rapgap were compared 
with other MC models: CASCADE; Rapgap with the Colour Dipole Model [40]; and 
Rapgap with the Peterson fragmentation function [41]. Two different values of the e 
parameter of the Peterson fragmentation function were used, namely e = 0.0055 and 
0.0041 as recently determined in e^e~ collisions by the SLD and OPAL collabora- 
tions, respectively [42]. The corresponding systematic uncertainty was defined as the 
maximal deviation with respect to the reference sample and was +2.2%. 

These systematic uncertainties were added in quadrature separately for the positive and 
negative variations to determine the overall systematic uncertainty. These estimates were 
also made in each bin in which the differential cross sections were measured. The uncer- 
tainty associated with the luminosity measurement for the 1999-2000 data-taking periods 
used in this analysis was ±2.2%. This introduces an overall normalisation uncertainty on 
each measured cross section, which is correlated between all data points. This is added 
in quadrature to the other systematic uncertainties on the total visible cross section, but 
is not included in the figures or tables of the differential cross section measurements. 

7 Results 

The total visible cross section, cr^^, was determined in the kinematic range > 2 GeV^, 
0.05 < y < 0.7 with at least one hadron-level jet in the Breit frame with E^^ct > 6 GeV 
and —2 < ?7jj^J' < 2.5 and with a muon fulfilling the following conditions: —0.9 < 77'* < 1.3 
and > 2 GeV or —1.6 < t]^ < —0.9 and > 2 GeV. The jets were defined by applying 
the /ct algorithm to stable hadrons; weakly decaying B hadrons are considered unstable. 
The muons coming from direct and indirect b decays are matched to any jet in the event. 
The measured cross section is 

^bbi^P — ^e66X— >-ejet/i X) = 40.9 ± 5.7 (stat.)^^ °(syst.) pb. 

This measurement has been corrected for electroweak radiative effects using Heracles. 
The NLO QCD prediction with hadronisation corrections is 20.6^2 2 pb, which is about 
2.5 standard deviations lower than the measured total cross section. The Cascade MC 
program gives cx^g = 28 pb and Rapgap gives a^g = 14 pb. 
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The differential cross sections were calculated in the same restricted kinematic range as 
the total cross section by repeating the fit of the p^rf^ distribution and evaluating the 
electroweak radiative corrections in each bin. The results are summarised in Table ^ 

FiguresEk andEb show the differential cross sections as functions of and x, respectively, 
compared to the NLO QCD calculation. The NLO QCD predictions generally agree with 
the data; in the lowest and lowest x bins, the data are about two standard deviations 
higher. Figures |2t and show the same differential cross sections compared with the 
Rapgap and Cascade MC simulations. Cascade agrees with the data except for the 
lowest and lowest x bin. Rapgap is well below the data in all bins, but it reproduces 
the shapes of the data distributions. 

Figures and Ob show the differential cross sections as functions of the transverse mo- 
mentum, p^, and pseudorapidity, rj^^, of the muon, compared to the NLO QCD calculation. 
They generally agree with the data; in the lowest bin and the high 77^ bin, the NLO 
QCD prediction is about two standard deviations below the data. FiguresEfc andEli show 
the same differential distribution compared with Cascade and Rapgap. Cascade de- 
scribes the measured cross sections well except for the lowest pj, bin, while Rapgap lies 
below the data. 

Figure EK shows the differential cross section as a function of -E^p jet* oi the leading jet 
compared to the NLO QCD calculation. The NLO QCD prediction agrees with the data 
reasonably well, though it is systematically below. For the highest E^Yct bin the difference 
is about two standard deviations. Figure shows the same differential distribution 
compared with Cascade and Rapgap. For all -E-^jct* values. Cascade reproduces the 
measured cross section reasonably well while Rapgap lies below the data. 



8 Conclusions 

The production of b quarks in the deep inelastic scattering process — > e /i jet X has 
been measured with the ZEUS detector at HERA. The NLO QCD prediction for the 
visible cross section lies about 2.5 standard deviations below the measured value. 

Single differential cross sections as functions of the photon virtuality, Q^, the Bjorken 
scaling variable, x, the transverse momentum and pseudorapidity of the muon as well 
as the transverse energy of the leading jet in the Breit frame have been measured. The 
Cascade MC program, implementing the CCFM QCD evolution equations, gives a good 
description of the measured cross sections. It is, however, below the data for low values 
of the transverse momenta, low and low values of x. Rapgap is well below the data 
for all measured cross sections. The differential cross sections are in general consistent 
with the NLO QCD predictions; however at low values of Q^, Bjorken x, and muon 
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transverse momentum, and high values of jet transverse energy and muon pseudorapidity, 
the prediction is about two standard deviations below the data 

In summary, 6-quark production in DIS has been measured for the first time and has been 
shown to be consistent with NLO QCD calculations. 
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range 
( GeV^) 


(pb/ GeV^) A,t,t ^,y,t 


2,10 
10,40 
40, 1000 


2.63 ±0.56 +!^;^^ 
0.36 ±0.10 

0.010 ±0.002 


log]^Q(x) range 


da jdx 

(pb) l^stat ^syst 


-4.5,-3.5 
-3.5,-2.9 
-2.9,-1.0 


20.9 ±4.4 tl'l 
17.2 ±4.7 tli 
5.3 ±1.3 1?:^ 


•p^ range 
(GeV) 


(icT / c/p^ 
(pb/ GeV) A^toj A.yst 


2,3 
3,4 

4,15 


30.5 ±7.6 tli 
9.7 ±2.6 111 
0.59 ±0.13 toll 


ri^ range 


da 1 di]^ 

(pb) l\^tat ^syst 


-1.6,-0.15 
-0.15,0.45 
0.45,1.3 


9.1 ±2.2 111 
14.2 ±3.6 tl'^Q 
19.8 ±4.1 tli 


E'^^^^l range 
(GeV) 


da/dE^^^ 

(pb/ GeV) Astat ^syst 


6,10 

10,13 

13,36 


5.7 ±1.4 +li 
3.4 ±0.8 
0.40 ±0.08 +°|^ 



Table 1: Single differential b- quark cross sections as functions of Q"^ , the Bjorken- 
X variable, the muon transverse momentum, pt^, the muon pseudorapidity, rj^, and 
the transverse energy of the leading jet in the Breit frame, E^^et- statistical 
and systematic uncertainties are shown separately. 
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Figure 1: (a) p^r^^ distribution measured for unidentified tracks in an inclusive 
DIS sample compared with the Rapgap MC simulation (see text). Data (dots) 
and the Rapgap MC (solid line) distributions after the final event selection for: 
(b) the measured p'^r^^ distribution; (c) muon momentum; (d) muon pseudorapidity; 
(e) transverse energy in the Breit frame; and (f) pseudorapidity in the laboratory 
frame of the associated jet. The solid line represents all MC contributions while the 
hatched histograms show the contribution from b quarks according to the percentage 
given by the fit (see Section^. The error bars are statistical only. 
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Figure 2: Differential b-quark cross section as a function of (a) and (b) 
Bjorken x for events with at least one jet reconstructed in the Breit frame and a 
muon, compared to the NLO QCD calculations. The error bars on the data points 
correspond to the statistical uncertainty (inner error bars) and to the statistical 
and systematic uncertainties added in quadrature (outer error bars). The solid 
line shows the NLO QCD calculations with the hadronisation corrections and the 
dashed line the same calculation without the hadronisation corrections. The shaded 
bands show the uncertainty of the NLO QCD prediction due to the variation of the 
renormalisation and factorisation scale, fi, and the b-quark mass, mi,. Differential 
b-quark cross sections as a function of (c) and (d) Bjorken x, compared with 
the LO QCD MC programs Cascade (solid line) and Rapgap (dashed line). 
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Figure 3: Differential b-quark cross section as a function of (a) the muon trans- 
verse momentum p'^ and (h) muon pseudorapidity rj^ in the laboratory frame, com- 
pared to the NLO QCD calculations. Other details are as described in the caption 
to Fig. Differential b-quark cross section as a function of (c) pi^ and (d) rj^ , 
compared with LO QCD MC programs CASCADE (solid line) and Rapgap (dashed 
line). 
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Figure 4: (a) differential b-quark cross section as a function of the transverse 
energy of the jet in the Breit frame E^^et- ^^^^ (dots) are compared to the 

NLO QCD calculations (a). Other details are as described in the caption to Fig.{^ 
(b) differential b-quark cross sections as a function of E^Yct compared with LO 
QCD MC programs Cascade (solid line) and Rapgap (dashed line). 
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